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Multinucleon transfer processes in heavy-ion reactions at energies slightly above the Coulomb bar- 
rier are investigated in a fully microscopic framework of the time-dependent Hartree-Fock (TDHF) 
theory. Transfer probabilities are calculated from the TDHF wave function after collision using the 
projection operator method which has recently been proposed by Simenel (C. Simenel, Phys. Rev. 
Lett. 105, 192701 (2010)). We show results of the TDHF calculations for transfer cross sections 
of the reactions of 40 ' 48 Ca+ 124 Sn at £ lab = 170, 174 MeV, 40 Ca+ 208 Pb at £ lab = 235, 249 MeV, 
and 58 Ni+ 208 Pb at _Ei a b = 328.4 MeV, for which measurements are available. We find the transfer 
processes show different behaviors depending on the N/Z ratios of the projectile and the target, 
and the product of the charge numbers, ZpZ^. When the projectile and the target have different 
N/Z ratios, fast transfer processes of a few nucleons towards the charge equilibrium of the initial 
system occur in reactions at large impact parameters. As the impact parameter decreases, a neck 
formation is responsible for the transfer. A number of nucleons are transferred by the neck breaking 
when two nuclei dissociate, leading to transfers of protons and neutrons in the same direction. Com- 
paring cross sections by theory and measurements, we find the TDHF theory describes the transfer 
cross sections of a few nucleons reasonably. As the number of transferred nucleons increases, the 
agreement becomes less accurate. The TDHF calculation overestimates transfer cross sections ac- 
companying a large number of neutrons when more than one proton are transferred. Comparing 
our results with those by other theories, we find the TDHF calculations give qualitatively similar 
results to those of direct reaction models such as GRAZING and Complex WKB. 



I. INTRODUCTION 

In the last three decades, measurements of multin- 
ucleon transfer process have been achieved extensively 
in heavy-ion collisions close to the Coulomb barrier [J- 
[24j ]. This process is an interesting quantum dynamics of 
many nucleons regarded as a nonequilibrium transport 
phenomenon which reflects both static properties and 
time-dependent dynamics of colliding nuclei. The static 
properties include nuclear shell structure and binding 
energies, nuclear shapes, properties of outermost single- 
particle orbitals, neutron-to-proton mass ratios, pairing 
properties, and so on. The dynamical effects include 
quantum tunneling between colliding nuclei, matching of 
Q-values and momenta, formation and breaking of the 
neck, quasi-fission dynamics, and so on. 

Besides fundamental interests in its mechanisms, the 
multinucleon transfer reaction is also expected to be use- 
ful as a mean to produce unstable nuclei whose produc- 
tion is difficult by other methods. For example, a pro- 
duction of neutron-rich nuclei with atomic mass around 
A ~ 200 along the neutron magic number N = 126 
has been discussed [25T - [28| . The knowledge on structural 
properties of these nuclei is crucially important to un- 
derstand a detail scenario of heavy elements synthesis in 
the r-process. An experiment to produce such neutron- 
rich unstable nuclei has been planned in the reactions of 
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Xe isotopes on 198 Pt (29J. A production of superheavy 
elements using multinucleon transfer reactions has also 
been discussed (see, for example, Refs. [28l l3fj| - [34j |). 

To describe multinucleon transfer processes theoret- 
ically, models based on a direct reaction picture such 
as GRAZING @ and Complex WKB (CWKB) M 
have been extensively developed and applied [H, [HI Il3h - 
[201 ] . In these models, multinucleon transfer processes are 
treated statistically, using single-nucleon transfer proba- 
bilities calculated by the first-order perturbation theory. 
A model based on Langevin-type equations of motion has 
also been developed [111 [HJ ■ This model describes not 
only multinucleon transfer processes but also deep in- 
elastic collisions, quasi-fission, fusion-fission, and fusion 
reactions in a unified way [13, [H, l3bTj33l | . 

Although the above mentioned approaches have shown 
reasonable successes, these models are not fully micro- 
scopic but include some model assumptions. To get fun- 
damental understanding of the dynamics and to present 
a reliable prediction for the cross sections, it is highly 
desired to develop a fully microscopic description for the 
multinucleon transfer processes with minimum assump- 
tions on the dynamics. To this end, we conduct micro- 
scopic calculations in the time-dependent Hartree-Fock 
(TDHF) theory. 

The theory of the TDHF was first proposed by Dirac 
in 1930 [37|. Applications of the TDHF theory to nuclear 
collision dynamics started in 1970s [38-45]. Progresses in 
the early stage have been summarized in Ref. [461 ] . Since 
then, continuous efforts have been devoted for improv- 
ing the method and extending applications [47H7l| . At 
present, three-dimensional calculations with full Skyrme 
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functionals including time-odd components are routinely 
conducted. In most TDHF calculations, Skyrme-type in- 
teractions [72] are used. Since parameters of Skyrme in- 
teractions are determined to reproduce nuclear properties 
for a wide mass region, there is no empirical parameter 
specific to the reaction. 

The TDHF theory may describe both peripheral and 
central collisions. In peripheral collisions, the mean-field 
of the collision partner works as a time-dependent per- 
turbation for the orbitals. This picture of the transfer dy- 
namics is similar to that in direct reaction models where 
single-particle transfer probabilities are calculated either 
by the perturbation theory [H, HH or by solving numer- 
ically the time-dependent Schrodinger equation (73l - [75| . 
In collisions at smaller impact parameters, the TDHF 
theory describes macroscopic dynamics such as fusion 
[H |H S3 M, M, H, HE Hf , ^asi-fissionJpj6lj65|, [n| , 
and deep inelastic collisions [34], [39M42L H . 14a . [6a | . Nu- 
cleons are exchanged between projectile and target nuclei 
through the neck formation. This description of the mult- 
inucleon transfer processes is similar to the Langevin- 
type description [3(3, [HJ ■ In this way, the TDHF theory 
is expected to be capable of describing quite different 
transfer mechanisms in a unified way. 

In this paper, we will apply the TDHF theory to calcu- 
late transfer probabilities and cross sections for the reac- 
tions of 40 Ca+ 124 Sn at E lab = 170 MeV, 48 Ca+ 124 Sn at 
E uh = 174 MeV, 40 Ca+ 208 Pb at E lah = 235, 249 MeV, 
and 58 Ni+ 208 Pb at E lah = 328.4 MeV, for which measure- 
ments are available [ltj 0, U3, Gil • To calculate transfer 
probabilities from the TDHF wave function after colli- 
sion, we use the projection operator method which has 
recently been proposed by Simenel [67j . 

In addition to the fact that extensive measurements 
are available for these systems, analyses and comparisons 
of these systems are of much interest since reactions in 
these systems are expected to show qualitatively differ- 
ent features. While 48 Ca+ 124 Sn has almost the same 
neutron-to-proton ratio, N/Z, between the projectile and 
the target, other three systems have different N/Z ra- 
tios. We expect transfer processes towards the charge 
equilibrium take place in collisions with large N/Z asym- 
metry [H, [70, [73. Moreover, it is well known that the 
basic feature of the low-energy heavy-ion collisions de- 
pends much on the product of the charge numbers of the 
projectile and the target nuclei, ZpZ^. Fusion reactions 
beyond the critical value, ZpZ^ ~ 1600, are known to ac- 
company an extra- push energy [34], [73] • Analyses of the 
fusion-hindrance phenomena in the TDHF theory have 
been reported in Ref. [34[, showing that the extra-push 
energy in the TDHF calculation is in good agreement 
with that of the Swiatecki's extra- push model [7g. The 
four systems to be analyzed have different ZpZ^ values, 
1000 for 40 > 48 Ca+ 124 Sn, 1640 for 40 Ca+ 208 Pb, and 2296 
for 58 Ni+ 208 Pb. 

It has been considered that the success of the TDHF 
theory is limited to observables expressed as expectation 
values of one-body operators. Indeed, the particle num- 
ber fluctuation in deep inelastic collisions has been found 



to be substantially underestimated in the TDHF calcu- 
lations [13, Hg. Since transfer probabilities in the 
TDHF calculation may not be given as expectation val- 
ues of any one-body operators, it is not at all obvious 
whether the TDHF calculation provides a reasonable de- 
scription for multinucleon transfer processes. One of the 
main purposes of the present paper is to clarify usefulness 
and limitation of the TDHF calculation for the multin- 
ucleon transfer processes. We note that Simenel has re- 
cently presented a calculation using the Barian-Veneroni 
prescription [79j and concluded that the particle number 
fluctuation may not be affected much by the correlation 
effects beyond the TDHF theory for reactions which are 
not so much violent as deep inelastic collisions |80| . 

The construction of this paper is as follows. In Section 
im we describe a formalism to calculate transfer proba- 
bilities from the TDHF wave function after collision. We 
also describe our computational method. In Section IlIIl 
we present results of our TDHF calculations for four sys- 
tems and compare them with measurements. In Section 
IIV1 we compare our results with those by other theories. 
In Section [Vj a summary and a future prospect will be 
presented. 



II. FORMULATION 

A. Definition of transfer probabilities 

We consider a collision of two nuclei described by the 
TDHF theory. The projectile is composed of Np nucle- 
ons and the target is composed of Nt nucleons. The 
total number of nucleons is N = N P + N T . In the TDHF 
calculation, a time evolution of single-particle orbitals, 
4>i(r,a,t) (i = 1, • ■ • , N), is calculated where r and a 
denote the spatial and the spin coordinates, respectively. 
The total wave function is given by the Slater determi- 
nant composed of the orbitals: 

$(x lr -- ,x N ,t) = -j= det{cj)i(x j: t)}, (1) 

where a; is a set of the spatial and the spin coordinates, 
x = (r, a). For the moment, we will develop a formalism 
for a many-body system composed of identical fermions. 
An extension to the actual nuclei composed of two kinds 
of fermions, protons and neutrons, is simple and obvious. 

Before the collision, two nuclei are separated spatially. 
We divide the whole space into two, the projectile re- 
gion, Vp, and the target region, Vj,. After the collision, 
we assume that there appear two nuclei, a projectile-like 
fragment (PLF) and a target-like fragment (TLF). We 
ignore channels in which nuclei are separated into more 
than two fragments after the collision. We again intro- 
duce a division of the whole space into two, the projectile 
region, Vp , which includes the PLF, and the target re- 
gion, , which includes the TLF. 

We define the number operator of each spatial region 
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as 



J T ■ i i 

l — l 1=1 

where r specifies the spatial region either Vp U) or V^ f) . 
We introduce the space division function, T (r), defined 
as 



6 T (r) = 



1 for r G r 
for r ^ t 



(3) 



The sum of the two operators, iV \ms) and N aj) , is the 

number operator of the whole space, N = N v ^ + N V i 

= N v f +N v f. In ordinary TDHF calculations, an initial 
p . T 

wave function is the direct product of the ground state 
wave functions of two nuclei boosted with the relative 
velocity. The single-particle orbitals, (f>i(x,t), are local- 
ized in one of the spatial regions, Vp or V^, at the initial 
stage of the calculation. Therefore, the initial wave func- 
tion is the eigenstate of both operators, N V i and Nyt, 
with eigenvalues, iVp and Nt, respectively. At the final 
stage of the calculation after the collision, each single- 
particle orbital extends spatially to both spatial regions 



of Vp and . Due to this fact, the Slater determinant 
at the final stage is not an eigenstate of the number oper- 
ators, N v f and N v f, but a superposition of states with 
different particle number distributions. 

The probability that n nucleons are in the spatial re- 
gion Vp and N — n nucleons are in the spatial region Vp 
is defined as follows. We start with the normalization 
relation of the final wave function after the collision, 
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dx N \^{xi, ■ ■ ■ ,x N )\ = 1, 
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where J dx = ^2 a J dr. Here and hereafter, we denote 
the many-body wave function at the final stage of the 
calculation as ty(xi,--- , x N ) = &et{ifji(xj)} /y/W\, and 
omit the time index. We also omit the suffix / from Vp 
and Vf. . The normalization relation, Eq. Q , includes N- 
fold integral over the whole spatial region. We divide each 
spatial integral into two integrals over the subspaces, Vp 
and Vt- We then classify the 2 N terms, generated by the 
divisions of the spatial regions, according to the number 
of Vp and the number of Vr included in the integral: 
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JVp+Vt 
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dx N \y(xi,- ■ ■ ,x N )\' 

Vp+V-r 



X! X! jdx v ■ ■jdx N \^{x l , 
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where each subscript r» (i = 1, ■ ■ • , N) represents either 
Vp or Vt. The notation s({n} ■ V^~ n ) means that 
the sum should be taken for all possible combinations 
of Ti on condition that, in the sequence of t\,--- ,tn, 



Vp appears n times and Vr appears N — n times. The 
number of the combinations equals to nC u . From this 
expression, we find the probability that n nucleons are in 
the Vp and N — n nucleons are in the Vt is given by 



P n = ^2 / dxi ■ ■ ■ I dx N \^(xi, 

, •> Tl J Tn 



■■■ ,X N ) 



(6) 
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Equation (|S|) ensures the relation, ^2 n=0 Pn — 1- From 
the probability P n , we may obtain nucleon transfer prob- 
abilities. For example, the probability of n-particle trans- 
fer from the projectile to the target is given by P/v P - n - 



B. Number projection operator 

Above expression of the probability P n can be repre- 
sented as an expectation value of the number projection 
operator P n , i.e. P n = (^\P n l^)- This operator extracts 
a component of the wave function with particle number 
n in the Vp and TV — n in the Vt from the final wave 
function ty(xx,--- ,xn). From Eq. ([5]), we obtain the 
following expression for the number projection operator, 

Pn= £ e T1 (ri)---e>T N (rN). (7) 

The projected wave function, P n ^, is the eigenstate of 
the number operators, N\/ P and Ny T , with eigenvalues, n 
and N — n, respectively. From Eq. ([3]), there follows 
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E^=n( e ^)+ e ^(n) =i- (8) 



Recently, Simenel has provided an alternative expres- 
sion for the number projection operator (oTj which is 
given by 
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Pn = 7T- / dBS n ~ Nv e 



(9) 



We can easily show that this expression, Eq. (j9|), is equiv- 
alent to Eq. ([7]) as follows: 
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C. Computation of transfer probabilities 

Two expressions for the number projection operator 
P n , Eq. (0 and Eq. (j9|), have been utilized to calculate 
transfer probabilities in the TDHF theory. When we use 
Eq. 0, the probability P n is expressed in terms of the 
single-particle orbitals as 



dxt ■ ■ J dx N ij}*{x\) ■ ■ ■ ^* N {x N )P n det{ipi(xj)} 



= J2 det{<^to) Tl }, 

s({r,}:V p "<-") 



(10) 



where the summation over £ is taken for all possible per- 
mutations of the index £j (i = 1, • • ■ , iV), and sgn(£) 
is a sign depending on the number of permutations. 
{ipi\ipj) T = J dx (x)ipj (x) denotes an overlap integral 
in the spatial region r. 

When we use Eq. ©, we obtain 
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Two expressions, Eq. (JTUJ) and Eq. ([IT]) , should give 
equivalent results. We indeed confirmed that both ex- 
pressions give the same results for light systems. How- 
ever, the computational cost is rather different between 
two methods. Let us first consider the computational 
cost of Eq. ([ID]). 

P n = det{(^) r J. 

s({T i }:V p »V*-") 

In this expression, it is necessary to calculate the determi- 
nants of dimension N many times. For example, to calcu- 
late the probabilities of all possible processes, Po to Pn, 
we need to calculate determinants of dimension N for 
2 N times. Even for the calculation of the probability 
without any particle transfer, we need to calculate the 
determinants as many as nCn p - The calculation in this 
way soon becomes impossible as N increases and is use- 
ful only for light systems. This method has been used 
in the 40 Ca+ 40 Ca collision in Ref. ■ It has also been 
used in the electron transfer processes in atomic collisions 

Earn. 

When we use the expression of Eq. , 
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the computational cost can be significantly small. In 
this expression, we achieve integral over 6 employing the 



trapezoidal rule discretizing the interval [0, 2ii\ into M 
equal grids. To calculate all the probabilities, Po to Pn, 
we need to calculate the determinants of dimension N 
for M times. We find M = 200 is sufficient for systems 
presented in this paper. In our calculations shown below, 
we employ Eq. (1X11) . 



D. Transfer cross sections 

We next derive the formula for cross sections of transfer 
reactions. We assume that both projectile and target 
nuclei are spherical, so that the reaction is specified by 
the incident energy E and the impact parameter b. 

Up to this point, we derived expressions of transfer 
probabilities for a system composed of identical fermions. 
Since the TDHF wave function is a direct product of 
Slater determinants for protons and neutrons, the reac- 
tion probability is also given by the product of the proba- 
bilities for protons and neutrons. Let us denote the prob- 
ability that Z protons are included in the Vp as P^ (b) 

and N neutrons are included in the Vp as P^' (b) . Then, 
the probability that Z protons and N neutrons are in- 
cluded in the Vp is given by 



Pz, N (b) = P < i\b)P { N\b). 



(12) 



We calculate the transfer cross section for the channel 
where the PLF is composed of Z, N nucleons by integrat- 
ing the probability Pz,n(P) over the impact parameter, 



a tI (Z,N) =2ir / bP ZiN (b)db 



(13) 



The minimum of the integration over the impact param- 
eter is the border dividing fusion and binary reactions. 
In practice, we first examine the maximum impact pa- 
rameter in which fusion reactions take place for a given 
incident energy. We will call it the fusion critical im- 
pact parameter and denote it as b{. We then repeat re- 
action calculations at various impact parameters for the 
region, b > bf, and calculate the cross section by numer- 
ical quadrature according to Eq. (fT3"|) . 



E. Numerical methods 

We have developed our own computational code of the 
TDHF theory for heavy-ion collisions. We employ a uni- 
form spatial grid in the three-dimensional Cartesian co- 
ordinate to represent single-particle orbitals without any 
symmetry restrictions. The grid spacing is taken to be 
0.8 fm. We take a box size of 60 x 60 x 26 grid points 
(48 fm x 48 fm x 20.8 fm) for collision calculations, 
where the reaction plane is taken to be the xy-plane. The 
initial wave functions of projectile and target nuclei are 
prepared in a box with 26 x 26 x 26 grid points. We use 
11-points finite-difference formula for the first and second 
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derivatives. To calculate the time evolution of single- 
particle orbitals, we use the Taylor expansion method 
of 4th order. The first-order predictor-corrector step is 
adopted in the time evolution. The time step is set to 
At = 0.2 fm/c. To calculate the Coulomb potential, we 
employ the Hockney's method [83[ in which the Fourier 
transformation is achieved in the grid of two times larger 
box than that utilized to express single-particle orbitals. 

We have tested the accuracy of the code by comparing 
our results with those by other codes. We have confirmed 
that the fusion critical impact parameters of the reac- 
tions of 16 0+ 16 and 16 0+ 28 reported in Ref. 01 are 
reproduced within 0.1 fm accuracy by our code. We have 
also calculated the fluctuation of exchanged nucleons for 
40 Ca+ 40 Ca head-on collisions and confirmed that results 
reported in Ref. [84| are reproduced accurately. 



III. RESULTS 

In this section, we will show calculated results for 
the reactions of 40 Ca+ 124 Sn at the incident energy of 
170 MeV, 48 Ca+ 124 Sn at 174 MeV, 40 Ca+ 208 Pb at 235, 
249 MeV, and 58 Ni+ 208 Pb at 328.4 MeV. 

As for the energy density functional and potential, we 
use the Skyrme functional including all time-odd terms 
except for the second derivative of the spin densities, 
As( n - p K We encounter numerical instability in the time 
evolution calculation if we include the term in the poten- 
tial. All of the results reported here are calculated using 
the Skyrme SLy5 parameter set [85| . This interaction 
has been utilized in the fully three-dimensional TDHF 
calculations for heavy- ion collisions |54j, |6l|, |69j . 

In the ground state calculations, we find the ground 
states of 40,48 Ca, and 208 Pb are spherical, the ground 
state of 124 Sn is almost spherical {(3 ~ 0.02), and the 
ground state of 58 Ni is slightly deformed (j3 ~ 0.1) in 
prolate shape. 

We take the incident direction parallel to the x-axis 
and the impact parameter vector parallel to the y-axis. 
The reaction is specified by the incident energy and the 
impact parameter. As an initial condition, the two col- 
liding nuclei are placed with the distance 16-18 fm in 
the x-direction. Before starting the TDHF calculation, 
we assume the centers of the two colliding nuclei follow 
the Rutherford trajectory. For the deformed nuclei, we 
placed the nucleus with the symmetry axis being set par- 
allel to the z-axis. 

We stop time evolution calculations when two nuclei 
are separated by 20-26 fm, if binary fragments are pro- 
duced. If the colliding nuclei fuse and do not separate, 
we continue time evolution calculations up to 1500 fm/c 
after two nuclei touch for 40 ' 48 Ca+ 124 Sn reactions, and 
more than 3000 fm/c after two nuclei touch for other 
reactions. 

For each collision system, we first find the fusion crit- 
ical impact parameter b{. We find them by repeating 
calculations changing the impact parameter by 0.01 fm 



step. We then calculate reactions for various impact pa- 
rameters outside the critical value. At an impact param- 
eter region smaller than 7 fm, we calculate reactions of 
impact parameters with 0.25 fm step. At an impact pa- 
rameter region larger than 7 fm, we calculate reactions 
of b = 7.5, 8, 9, and 10 fm. Close to the fusion critical 
impact parameter, we calculate reactions in 0.05 fm and 
0.01 fm impact parameter steps. All these calculations 
are used to evaluate the transfer cross sections. In calcu- 
lating transfer cross sections according to Eq. (fT5")) . the 
upper limit of the integral over b is set to 10 fm. 

A. 40 ' 48 Ca+ 124 Sn reactions 

In this subsection, we present results for the reac- 
tions of 40 Ca+ 124 Sn at E l&h = 170 MeV (£ c . m . ~ 
128.5 MeV) and 48 Ca+ 124 Sn at E lah = 174 MeV (E c . m . ~ 
125.4 MeV), for which multinucleon transfer cross sec- 
tions have been measured experimentally [l(| EH . The 
neutron-to-proton ratio, N/Z, is different between the 
projectile and the target for 40 Ca+ 124 Sn, while it is al- 
most the same for 48 Ca+ 124 Sn. Therefore, we may ex- 
pect different features in the transfer process. As we men- 
tioned in the introduction, the product of charge numbers 
of the projectile and the target is important for the fu- 
sion dynamics. The present systems have ZpZ^ = 1000 
< 1600, so that no fusion-hindrance is expected to occur. 

To estimate the Coulomb barrier height, we cal- 
culate the nucleus-nucleus potential using the frozen- 
density approximation neglecting the Pauli blocking ef- 
fect [63, il, H3j. The potential is given by V(R) = 
E[p r + p T ](R) — -Eg.s. [p P ] — -Eg.s. [p T ], where R is the dis- 
tance between the centers-of-masses of the two nuclei, 
p p (p T ) denotes nuclear density of the projectile (target) 
in their ground state. E[p p + p T ](R) denotes the to- 
tal energy when two nuclei are separated by the relative 
distance R. -E g . s . [p P ] and -E g . s . [p T ] denote the ground 
state energy of each nucleus. In the calculation, the 
Coulomb barrier height is estimated as Vb ~ 118.2 MeV 
for 40 Ca+ 124 Sn and V B « 115.8 MeV for 48 Ca+ 124 Sn, 
respectively. Since the initial relative energies are higher 
than the Coulomb barrier heights, we find the fusion crit- 
ical impact parameter, b{ = 3.69 fm for 40 Ca+ 124 Sn and 
6f = 3.62 fm for 48 Ca+ 124 Sn, respectively. 

1. Overview of the reactions 

Before showing detailed analyses of transfer reactions, 
we first present an overview of the reaction dynamics. 
In Fig. Q] we show the deflection function, 9(6), in (a) 
and the total kinetic energy loss (TKEL) in (b), as func- 
tions of impact parameter b. Results for the 40 Ca+ 124 Sn 
reactions are denoted by red filled triangles connected 
with solid lines, while results for the 48 Ca+ 124 Sn reac- 
tions are denoted by green open circles connected with 
dotted lines. In Fig. [T] (a), we also show deflection func- 
tions of the pure Coulomb trajectories by a red dotted 
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FIG. 1. (Color online) Deflection function (a) and total ki- 
netic energy loss (b) as functions of impact parameter 6 for the 
reactions of 40 Ca+ 124 Sn at £i ab = 170 MeV and 48 Ca+ 124 Sn 
at Ei ah = 174 MeV. Results for the 40 Ca+ 124 Sn reactions 
are denoted by red filled triangles connected with solid lines, 
while results for the 48 Ca+ 124 Sn reactions are denoted by 
green open circles connected with dashed lines. In (a), we 
also show deflection functions for the pure Coulomb trajecto- 
ries by a red dotted line for the 40 Ca+ 124 Sn reactions and by 
a green two-dot chain line for the 48 Ca+ 124 Sn reactions. 



line for 40 Ca+ 124 Sn and by a green two-dot chain line for 
48 Ca+ 124 Sn. 

The TKEL is evaluated from the center-of-mass mo- 
tions of the two fragments at the final stage of the cal- 
culation. Denoting the mass, charge number, and the 
momentum of the PLF (TLF) at the final stage of the 
calculation as M PLF ( TLF ), Z PL f(tlf), and P PL f(tlf), 
the TKEL is evaluated by TKEL = E c . m . - [(Pplf - 
-Ptlf) 2 /2/i + Z PLF ZTL F e 2 /Rf] . Here E c . m . is the initial 
incident energy in the center-of-mass frame, fi is the re- 
duced mass given by /z = Mplf-Mtlf/C^plf + -^tlf), 
and Rf is the distance between the centers-of- masses of 
the PLF and the TLF. The TKEL increases rapidly as 
the impact parameter decreases in the region b < 4.5 fm, 
where the deflection function, 0(6), decreases apprecia- 
bly by the nuclear attractive interaction. 

The deflection function shows a maximum at b ~ 4 fm 
and decreases inside this impact parameter. The maxi- 
mum deflection angle corresponds to the Coulomb rain- 
bow angle, 9 T . It is given by 101° for 40 Ca+ 124 Sn and 
103° for 48 Ca+ 124 Sn. In Fig.0 we compare the Coulomb 
rainbow angle for the 40 Ca+ 124 Sn reaction with mea- 
sured differential cross sections reported in Ref . [l(| ■ Red 
filled circles denote measured cross sections and blue solid 
vertical lines denote the Coulomb rainbow angle in the 
laboratory frame. As seen from the figure, the peak posi- 
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FIG. 2. (Color online) Differential cross sections of repre- 
sentative transfer channels as functions of scattering angle in 
the laboratory frame for the 40 Ca+ 124 Sn reaction at £"i a b = 
170 MeV. The Coulomb rainbow angle obtained from the 
TDHF trajectories is denoted by blue solid vertical lines, 
and is compared with measured differential cross sections, red 
filled circles, which have been reported in Ref. [lOj. 



tions of the measured cross sections roughly coincide with 
the Coulomb rainbow angle by the TDHF calculation. 

In Fig. [3J we show snapshots of density distribution for 
the 40 Ca+ 124 Sn reaction at b = 3.70 fm, just outside the 
fusion critical impact parameter, bf . We find a formation 
of a neck between the projectile and the target during 
the collision. As will be shown later, several nucleons are 
exchanged between the projectile and the target at this 
impact parameter. We find a formation of the neck for 
the impact parameter region smaller than b ~ 4 fm where 
the TKEL becomes appreciable. 

We next consider the average number of transferred nu- 
cleons and its fluctuation. We denote the average number 
of nucleons in the PLF as (q = n for neutrons, p 

for protons), which is calculated from the density distri- 
bution at the final stage of the calculation, 



(?) 

PLF 



dr p (q) {r), 

ound PLF 



(14) 



where p^(r) is the density distribution of neutrons 
(q = n) or protons (q = p). The spatial integration is 
achieved over a sphere whose center coincides with the 
center-of-mass of the PLF. The radius of the sphere is 
taken to be 10 fm. We calculate the average number of 
nucleons in the TLF in the same way taking the radius 
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FIG. 3. (Color online) Snapshots of density distribution of 
the 40 Ca+ 124 Sn reaction at £a ab = 170 MeV and b = 3.70 fm, 
just outside the fusion critical impact parameter. 
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FIG. 4. (Color online) Left panels for 40 Ca+ 124 Sn at E lah = 170 MeV and right panels for 48 Ca+ 124 Sn at £a ab = 174 MeV. 
(a) and (e): Average number of transferred nucleons from the target to the projectile, (b) and (f): Neutron-to-proton ratios, 
N/Z, of the PLF and the TLF after collision, (c) and (g): Average number of nucleons emitted to the continuum, (d) and 
(h): Fluctuation of transferred nucleon number. The horizontal axis is the impact parameter b. In (b) and (f), the equilibrium 
N/Z value of the total system is indicated by a horizontal dashed line. 



of 14 fm for the TLF. We summarize various expressions 
for the average number and the fluctuation of transferred 
nucleons in Appendix lAl 

We denote the neutron (proton) number of the projec- 
tile and the target as Np^ and , respectively. In gen- 
eral, there holds N^ F + iV^ F < N p q) + N^ q) , since some 
nucleons are emitted to the continuum by the breakup 
process. As will be shown later, however, the number of 
nucleons emitted to the continuum is very small in the 
present calculations. The average number of transferred 
nucleons from the projectile to the target, iV t r , is given 



by 

<W=ArW F -ArW. (15) 

Figure 3] shows the average number of transferred nu- 
cleons, iv£ , m ( a ) an d ( e )j the neutron-to-proton ra- 
tios, N/Z, of the PLF and the TLF after collision in 
(b) and (f), the average number of nucleons emitted to 
the continuum in (c) and (g), and the fluctuation of the 
transferred nucleon number in (d) and (h), as functions 
of impact parameter b for 40 ' 48 Ca+ 124 Sn reactions. 

In Fig. 0] (a) and (e) , the average number of transferred 
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neutrons is shown by filled symbols connected with solid 
lines, while the average number of transferred protons is 
shown by open symbols connected with dotted lines. Pos- 
itive values indicate the increase of the projectile nucleons 
(transfer from 124 Sn to 40 ' 48 Ca) and negative values in- 
dicate the decrease (transfer from 40,48 Ca to 124 Sn). As 
seen from Fig.|4](a) and (e), a large value of average num- 
ber of transferred nucleons is seen for 40 Ca+ 124 Sn at the 
impact parameter region close to the fusion critical im- 
pact parameter, while the average number of transferred 
nucleons is small for 48 Ca+ 124 Sn. 

We show in Fig. [4] (b) and (f) the neutron-to-proton 
ratios, N/Z, of the PLF and the TLF. For the PLF, it is 
given by N$ F /N<g F , and for the TLF by N$ F /N$> F . 
Before the collision, the N/Z ratio is given by 1.00 for 
40 Ca, 1.40 for 48 Ca, and 1.48 for 124 Sn. The N/Z ratio 
of the PLF (TLF) is denoted by filled (open) symbols con- 
nected with solid (dotted) lines. We also denote the N/Z 
ratio of the total system by a horizontal dashed line in the 
figure, 1.34 for 40 Ca+ 124 Sn and 1.46 for 48 Ca+ 124 Sn. We 
find the nucleons are transferred towards the direction of 
the charge equilibrium. Namely, protons are transferred 
from 40 Ca to 124 Sn, while neutrons are transferred from 
124 Sn to 40 Ca in the 40 Ca+ 124 Sn reaction. The N/Z ra- 
tios of the projectile and the target do not differ much 
for 48 Ca+ 124 Sn, and we find a small number of trans- 
ferred nucleons on average for this reaction. The average 
number of transferred nucleons decreases rapidly as the 
impact parameter increases. For the impact parameter 
region larger than b ~ 6 fm, the average number of trans- 
ferred nucleons almost vanishes. 

In Fig. [4] (c) and (g), we show the average num- 



ber of nucleons emitted to the continuum, TV, 



(?) 

break-up 



Ng J ) - {N P Y F + Ntlf)' durin S the time evolu " 
tion. The average number of neutrons (protons) emitted 
to the continuum is denoted by filled (open) symbols con- 
nected with solid (dotted) lines. As seen in the figure, the 
number of emitted nucleons is very small. The maximum 
value, about 0.02, is seen at the impact parameter close 
to the fusion critical impact parameter. 

In Fig. @] (d) and (h), we show the fluctuation of the 
transferred nucleon number. The expression for the fluc- 
tuation is given by Eq. (IA4[) . The fluctuation of the 
transferred neutron (proton) number is denoted by filled 
(open) symbols connected with solid (dotted) lines. We 
find the fluctuation decreases as the impact parameter 
increases. The fluctuation decreases more slowly than 
the average number of transferred nucleons as a func- 
tion of impact parameter. We also find the fluctuation of 
48 Ca+ 124 Sn is somewhat smaller than but comparable in 
magnitude to that of 40 Ca+ 124 Sn, although the average 
number is vanishingly small for 48 Ca+ 124 Sn. 



2. Transfer probabilities 

We next show transfer probabilities as functions of 
impact parameter which are obtained from the final 
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FIG. 5. (Color online) Neutron (left panels) and proton (right 
panels) transfer probabilities as functions of impact parame- 
ter b for the reactions of 40 Ca+ 124 Sn at E lllb = 170 MeV 
(upper panels) and 48 Ca+ 124 Sn at E lah = 174 MeV (lower 
panels). The positive (negative) number of transferred nu- 
cleons represents the number of nucleons added to (removed 
from) the projectile. Shaded regions at small impact pa- 
rameter (6 < 3.69 fm for 40 Ca+ 124 Sn and b < 3.62 fm for 
48 Ca+ 124 Sn) correspond to the fusion reactions. 



wave functions using the number projection procedure of 
Eq. (ITTj) . The nucleon transfer probabilities, P« (6), are 
shown in Fig. [5] (linear scale) and in Fig. [5] (logarithmic 
scale). In these figures, shaded regions at small impact 
parameter (b < 3.69 fm for 40 Ca+ 124 Sn and b < 3.62 fm 
for 48 Ca+ 124 Sn) correspond to the fusion reactions. The 
positive (negative) number of transferred nucleons repre- 
sents the number of nucleons added to (removed from) 
the projectile. Left panels show transfer probabilities of 
neutrons, while right panels show transfer probabilities 
of protons. 

From the figure, we find that probabilities of single- 
nucleon transfer (green dashed lines) extend to a large 
impact parameter region. As the number of transferred 
nucleons increases, the reaction probability is sizable only 
at a small impact parameter region, close to the fusion 
critical impact parameter. 

The directions of the transfer processes are the same as 
those we observed in the average number of transferred 
nucleons in Fig. [3] (a) and (e). Namely, in the case of 
40 Ca+ 124 Sn (upper panels of Fig. [SJ, protons are trans- 
ferred from 40 Ca to 124 Sn and neutrons are transferred 
from 124 Sn to 40 Ca, the directions towards the charge 
equilibrium. We note that the transfer probabilities to- 
wards the opposite directions, proton transfer from 12 Sn 
to 40 Ca and neutron transfer from 40 Ca to 124 Sn, are 
very small and are hardly seen in the linear scale figure 
(Fig. [5]). In the logarithmic scale (upper panels of Fig. [5]), 
we find the transfer probabilities towards the opposite di- 
rection to the charge equilibrium are smaller than those 
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FIG. 6. (Color online) Transfer probabilities in Fig. [5] are 
shown in logarithmic scale. Nucleon transfer probabilities op- 
posite to the direction of the charge equilibrium, which are 
not included in Fig. [5J are shown as well. 



towards the charge equilibrium by at least an order of 
magnitude. In the case of 48 Ca+ 124 Sn reaction (lower 
panels of Figs. [5] and [5]), the transfer probabilities towards 
both directions are the same order of magnitude. This 
is consistent with the fact that the average number of 
transferred nucleons is very small as shown in Fig. 0] (e). 



3. Transfer cross sections 

Integrating the transfer probabilities over impact pa- 
rameter, we obtain transfer cross sections. The results 
are shown in Fig. [7J for 40 Ca+ 124 Sn and in Fig. [5] for 
48 Ca+ 124 Sn. 

We first examine the 40 Ca+ 124 Sn reaction. Figure [7] 
shows the transfer cross sections classified according to 
the change of the proton number of the PLF from 40 Ca, 
as functions of neutron number of the PLF. Red filled 
circles denote measured cross sections and red solid lines 
denote results of the TDHF calculations. We show trans- 
fer cross sections of one proton added to (+lp) through 
six proton removed from (— 6p) 40 Ca. 

For zero- or one-proton removal channels, (Op) and 
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FIG. 7. (Color online) Cross sections for transfer channels 
classified according to the change of the proton number of the 
PLF from 40 Ca, as functions of neutron number of the PLF 
for the 40 Ca+ 124 Sn reaction at £i ab = 170 MeV. Red filled cir- 
cles denote measured cross sections and red solid lines denote 
results of the TDHF calculations. The number of transferred 
protons is indicated as (a;p) (—6 < x < +1). The measured 
cross sections have been reported in Ref. [lOI ] . 



(— lp), we find the experimental data are quantitatively 
reproduced by the TDHF calculations for up to five- 
neutron transfer processes. The calculated cross sections 
towards the direction opposite to the charge equilibrium 
are small, consistent with the observation in transfer 
probabilities shown in the upper panels of Fig. [51 

As the number of transferred protons increases, there 
appear some discrepancies between the TDHF calcula- 
tions and the measurements. The TDHF cross sections 
become too small compared with the measurements. We 
also find a shift of the peak position towards the larger 
neutron number. 

In Ref. (lOI . cross sections calculated by the GRAZ- 
ING code [35[ were compared with the measurements. 
In the GRAZING calculation, a similar discrepancy was 
observed. As the origin of the discrepancy, the signifi- 
cance of the evaporation effects has been mentioned [l(| ■ 
We will compare our results with those of the GRAZING 
calculations in Sec. IIVI 

We note that particle evaporation processes are not 
taken into account sufficiently in the present calculation. 
In Fig. [1] (b), we find the TKEL of as large as 25 MeV at 
a small impact parameter region where appreciable mult- 
inucleon transfer probabilities are found. The amount of 
the TKEL is sufficiently large to emit some nucleons to 
the continuum. However, as we saw in Fig. [4] (c), the 
average number of nucleons emitted to the continuum is 
very small, the maximum value is only 0.02. Although we 
have not yet estimated the number of evaporated nucle- 
ons, the inclusion of the evaporation processes is expected 
to reduce the discrepancy as follows. Neutron evapora- 
tion processes will shift the peak position of the transfer 
cross sections towards the smaller neutron number (left 
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FIG. 8. (Color online) Cross sections for transfer channels 
classified according to the change of the proton number of the 
PLF from 48 Ca, as functions of neutron number of the PLF 
for the 48 Ca+ 124 Sn reaction at £i ab = 174 MeV. Red filled cir- 
cles denote measured cross sections and red solid lines denote 
results of the TDHF calculations. The number of transferred 
protons is indicated as (xp) ( — 2<x< +2). The measured 
cross sections have been reported in Ref. [13| . 



direction in Fig. [7]). We may also expect that proton 
evaporation processes will shift cross section of n-proton 
removal channels to (n + l)-proton removal channels. 

Figure [5] shows transfer cross sections of 48 Ca+ 124 Sn 
reaction. The cross sections obtained from the TDHF 
calculations are in good agreement with the experimen- 
tal data for zero- and one-proton transfer channels, (Op) 
and (±lp). For two- proton transfer channels (±2p), how- 
ever, our TDHF calculations underestimate the cross sec- 
tions. In the case of two-proton removal channels (— 2p), 
the peak position shifts towards larger neutron number, 
while in the case of two-proton pickup channels (+2p), 
the peak position shifts towards smaller neutron num- 
ber. The underestimation in the (— 2p) channels may be 
remedied by taking into account the neutron evaporation 
processes as in the case of 40 Ca-|- 124 Sn reaction. However, 
the underestimation in the (+2p) channels may not. A 
similar discrepancy was reported in the GRAZING calcu- 
lation [l3j]. In Ref. [l3|], more complex mechanisms such 
as neutron- proton pair transfer and/or a-cluster transfer 
have been advocated for the origin of the discrepancy. 



B. 



) Ca+ 208p b reaction 



In this subsection, we present results for the reactions 
of 40 Ca+ 208 Pb at E lah = 235 and 249 MeV (E c . m . ~ 
197.1 and 208.8 MeV), for which measurements have been 
reported in Ref. [ljj]. This system has ZpZ^ = 1640, 
close to 1600. Therefore, we expect an appearance of 
the indication of the fusion-hindrance. We estimate the 
Coulomb barrier height of this system using the frozen- 
density approximation, giving Vb ~ 178.4 MeV. Since the 
collision energies are higher than the barrier height, we 
find finite values of the fusion critical impact parameter 
bf, as in the 40 ' 48 Ca+ 124 Sn reactions. They are given 
by h{ = 3.81 and 4.55 fm at £ lab = 235 and 249 MeV, 
respectively. 



1. Overview of the reactions 

We first present an overview of the reaction dynam- 
ics. In Fig. [SJ we show the deflection function in (a) and 
the TKEL in (b), as functions of impact parameter. Re- 
sults for the reaction at _Ei a b — 235 MeV are denoted 
by red filled triangles connected with solid lines, while 
results for the reaction at E\ a \-, = 249 MeV are denoted 
by green open circles connected with dotted lines. In 
Fig. [9] (a), we also show deflection functions for the pure 
Coulomb trajectories at £i a b = 235 MeV by a red dotted 
line and at -Ei a b = 249 MeV by a green two-dot chain 
line. In this system, we find an increase of the TKEL up 
to around 50 MeV and 60 MeV for the incident energies 
of 235 MeV and 249 MeV, respectively. This maximum 
value of TKEL is about a factor of two larger than the 
case of 40:48 Ca-f- 124 Sn reactions. We find the difference 
of the TKEL between these systems, 40 Ca+ 208 Pb and 
40 ' 48 Ca+ 124 Sn, comes from properties of the neck whose 
formation is observed when the TKEL becomes substan- 
tial. 

In Fig. [TUl we show snapshots of density distribution 
for the 40 Ca+ 208 Pb reaction at E uh = 249 MeV and b = 
4.56 fm, just outside the fusion critical impact parameter. 
The neck is seen to be formed solidly for a long period 



40 Ca + 208 Pb 




b(fm) 

FIG. 9. (Color online) Deflection function (a) and total ki- 
netic energy loss (b) as functions of impact parameter b for 
the reactions of 40 Ca+ 208 Pb at £i ab = 235 and 249 MeV. 
Results for the reactions at _Ei ab = 235 MeV are denoted by 
red filled triangles connected with solid lines, while results for 
the reactions at £7i a b = 249 MeV are denoted by green open 
circles connected with dashed lines. In (a), we also show de- 
flection functions for the pure Coulomb trajectories at _E] ab = 
235 MeV by a red dotted line and at _Bi ab = 249 MeV by a 
green two- dot chain line. 
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FIG. 10. (Color online) Snapshots of density distribution of 
the 40 Ca+ 208 Pb reaction at E iab = 249 MeV and b = 4.56 fm, 
just outside the fusion critical impact parameter. 

from t = 200 fm/c to 3000 fm/c. This process may be 
regarded as a quasi-fission. As will be shown below, a 
number of nucleons are transferred from 208 Pb to 40 Ca 
at this impact parameter. 

The period of the neck formation is longer in the 
present 40 Ca+ 208 Pb case than that in the 40 < 48 Ca+ 124 Sn 
cases. We find the neck formation for the periods of 1000- 
3000 fm/c and more for the present system depending on 
the impact parameter, while it is at most 300 fm/c in 
the 40,48 Ca+ 124 Sn systems. We consider this difference 
is related to the different ZpZ^ values of these systems. 
Since Z-pZ^ > 1600 in the present system, fusion reac- 
tions are hindered by the quasi-fission process. Namely, 
there appears a certain impact parameter region in which 
binary final fragments are produced after a rather solid 
neck formation during the collision. 

The Coulomb rainbow angle is 8 r ~ 99° for the reac- 
tion at £i a b = 235 MeV and 9 Y ~ 86° for the reaction 
at Ei a b = 249 MeV, respectively. The deflection func- 
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FIG. 11. (Color online) Differential cross sections of represen- 
tative transfer channels as functions of scattering angle in the 
center-of-mass frame for the 40 Ca+ 208 Pb reactions at _Bi a b = 
235 and 249 MeV. The Coulomb rainbow angle obtained from 
the TDHF trajectories is denoted by red solid (green dotted) 
vertical lines for E\ a b — 235 (249) MeV. They are compared 
with measured differential cross sections, red filled triangles 
(green open circles) for _E"i a b = 235 (249) MeV, which have 
been reported in Ref. 



tion becomes negative at small impact parameter region, 
reaching —200° at the fusion critical impact parameter. 
In Fig. 111! we compare the Coulomb rainbow angles for 
the 40 Ca+ 208 Pb reactions at E lah = 235 and 249 MeV 
with measured differential cross sections which have been 
reported in Ref. [l9|]. Red filled triangles denote mea- 
sured cross sections for E\ a b = 235 MeV, while green open 
circles denote those for E^ = 249 MeV. The Coulomb 
rainbow angle obtained from the TDHF trajectories is de- 
noted by red solid (green dotted) vertical lines for — 
235 (249) MeV. We find the peak positions of measured 
angular distributions are reasonably reproduced by the 
TDHF calculation. 

Figure [T2l shows the average number of transferred nu- 
cleons in (a), the N/Z ratios of the PLF and the TLF in 
(b), the average number of nucleons emitted to the con- 
tinuum in (c), and the fluctuation of the transferred nu- 
cleon number in (d), as functions of impact parameter. In 
each panel, triangles represent results for £q a b = 235 MeV 
and circles represent results for E\ a ^ — 249 MeV. 

In Fig. [T2] (a), the average number of transferred neu- 
trons is shown by filled symbols connected with solid 
lines, while the average number of transferred protons 
is shown by open symbols connected with dotted lines. 
Positive values indicate the increase of the projectile nu- 
cleons (transfer from 208 Pb to 40 Ca) and negative values 
indicate the decrease (transfer from 40 Ca to 208 Pb). As 
seen from the figure, the average number of transferred 
protons shows a minimum at a certain impact parame- 
ter {b = 4.0 fm for E lah = 235 MeV and b = 5.0 fm for 
Eiab = 249 MeV). Outside this impact parameter, the nu- 
cleon transfer process proceeds towards the direction of 
the charge equilibrium of the projectile and the target. 
Inside this impact parameter, neutrons are still trans- 
ferred towards the same direction. However, the number 
of transferred protons decreases and becomes positive, 
which corresponds to the transfer from 208 Pb to 40 Ca. 

At first sight, the direction of the proton transfer at 
small impact parameter region is opposite to the direc- 
tion of the charge equilibrium. However, it is not the case 
as can be understood from Fig. [T2] (b) which shows the 
neutron-to-proton ratios, N/Z, of the PLF (filled symbols 
connected with solid lines) and the TLF (open symbols 
connected with dotted lines), which are obtained from 
the average numbers of the nucleons shown in Fig. [T^] (a) . 
Before collision, the N/Z ratio is given by 1.00 for 40 Ca 
and 1.54 for 208 Pb. In Fig. [lj (b), the N/Z ratio of the 
total system, 1.43, is shown by a horizontal dashed line. 
As seen from the figure, the nucleon transfer processes 
proceed towards the direction of the charge equilibrium 
for both the PLF and the TLF at all impact parame- 
ter region outside the fusion critical impact parameter. 
Even though the average number of transferred protons 
shows complex behavior at small impact parameter re- 
gion, the N/Z ratios of the PLF and the TLF monoton- 
ically approach to the fully equilibrated value of 1.43 as 
the impact parameter decreases. 

The change of sign of the average number of transferred 
protons at small impact parameter region is found to be 
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Ca+ U0 Pb (£ lab =235, 249 MeV) 




ft (fin) 

FIG. 12. (Color online) The 40 Ca+ 208 Pb reactions at E lab = 
235 and 249 MeV. (a): Average number of transferred nucle- 
ons from the target to the projectile, (b): Neutron-to-proton 
ratios, N/Z, of the PLF and the TLF after collision, (c): 
Average number of nucleons emitted to the continuum, (d): 
Fluctuation of transferred nucleon number. The horizontal 
axis is the impact parameter b. Results for the reactions at 
-Blab = 235 MeV are denoted by triangles, while results for 
the reactions at Eiab = 249 MeV are denoted by circles. In 
(b), the equilibrium N/Z value of the total system, 1.43, is 
indicated by a horizontal dashed line. 

related to the formation of a rather solid neck. When 
the neck is broken, we find that most part of the neck 
is absorbed by the lighter fragment (c/. Fig. HU]). Since 
the neck is composed of both neutrons and protons, the 
absorption of the nucleons in the neck region results in 
the increase of average number of nucleons in the PLF 
for both neutrons and protons (see Fig. [T2l (a)). 



In Fig. [12] (c), we show the average number of nucle- 
ons emitted to the continuum during the time evolution. 
The average number of neutrons (protons) emitted to the 
continuum is denoted by filled (open) symbols connected 
with solid (dotted) lines. We count it by subtracting the 
number of nucleons inside a sphere of 14 fm for the TLF 
and that inside a sphere of 10 fm for the PLF from the 
total number of nucleons, 248. The average number of 
emitted nucleons is again very small, at most 0.1 around 
the fusion critical impact parameter. 

In Fig. [T^] (d), we show the fluctuation of the trans- 
ferred nucleon number. The fluctuation of transferred 
neutron (proton) number is denoted by filled (open) sym- 
bols connected with solid (dotted) lines. The fluctua- 
tion increases monotonically as the impact parameter de- 
creases, reaching the maximum value roughly 1.3 around 
the fusion critical impact parameter. Although the aver- 
age number of transferred protons is small at the small 
impact parameter region, the fluctuation of transferred 
proton number has value as large as that of neutrons. 
This fact indicates that single-particle wave functions 
of protons are exchanged actively between the projec- 
tile and the target, although the number of transferred 
protons is small on average. 



2. Transfer probabilities 

The nucleon transfer probabilities, Pn(b), are shown 
in Fig. [T31 In the figure, shaded regions at small im- 
pact parameter (6 < 3.81 fm for E\ a ^ = 235 MeV and 
b < 4.55 fm for £"i a b = 249 MeV) correspond to the fusion 
reactions. The positive (negative) number of transferred 
nucleons represents the number of nucleons added to (re- 
moved from) the projectile. In the left panels, we show 
transfer probabilities, (On) to (+4n) for neutrons and (Op) 
to (— 4p) for protons. In the right panels, we show trans- 
fer probabilities, (+ln) to (+9n) for neutrons and (— 4p) 
to (+4p) for protons, at small impact parameter regions 
just outside the fusion critical impact parameter. Prob- 
abilities of neutron transfer from 40 Ca to 208 Pb are very 
small and are not shown. 

As in the case of 40 ' 48 Ca-|- 124 Sn reactions, we find 
that probabilities of single-nucleon transfer (green dashed 
lines) extend to a large impact parameter region. Re- 
action probabilities for multinucleon transfer processes 
become appreciable at a small impact parameter region 
close to the fusion critical impact parameter. The trans- 
fer probabilities towards the charge equilibrium are large 
in most cases. At a small impact parameter region just 
outside the fusion critical impact parameter, however, we 
find substantial probabilities for the proton transfer pro- 
cesses opposite to the charge equilibrium as seen in the 
right panels of Fig. [T3] This is related to the increase 
of the average number of transferred protons at small 
impact parameter region which was seen in Fig. [T^] (a). 
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FIG. 13. (Color online) Neutron and proton transfer proba- 
bilities as functions of impact parameter b for the reactions of 
40 Ca+ 208 Pb at £ lab = 235 MeV (upper panels) and 249 MeV 
(lower panels). The positive (negative) number of transferred 
nucleons represents the number of nucleons added to (removed 
from) the projectile. Note that horizontal scales are differ- 
ent between the left and the right panels. Shaded regions at 
small impact parameter (b < 3.81 fm for = 235 MeV and 
b < 4.55 fm for E\ ah — 249 MeV) correspond to the fusion 
reactions. 



3. Transfer cross sections 

We show transfer cross sections in Figs. [14] and [15J 
Each panel of Fig. [H] shows cross sections classified ac- 
cording to the change of the proton number of the PLF 
from 40 Ca which is indicated by (xp) (—6 < x < +5), as 
functions of neutron number of the PLF. Each panel of 
Fig. [15] shows cross sections classified according to the 
change of the neutron number of the PLF from 40 Ca 
which is indicated by (xn) (—5 < x < +9), as functions 
of proton number of the PLF. Red filled triangles denote 
measured cross sections for E\ a i, = 235 MeV, while green 
open circles denote those for £"i a b = 249 MeV. Cross sec- 
tions calculated by the TDHF are denoted by red solid 
(green dotted) lines for £q ab = 235 (249) MeV. As seen in 
the average number of transferred nucleons in Fig. [T^] (a) 
and in the transfer probabilities in Fig. Q21 the transfer 
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FIG. 14. (Color online) Transfer cross sections for the 
40 Ca+ 208 Pb reactions at E ub = 235 and 249 MeV. Red filled 
triangles (green open circles) denote measured cross sections 
at £i a b = 235 (249) MeV. Red solid (green dotted) lines de- 
note results of the TDHF calculations at _Bi ab = 235 (249) 
MeV. The number of transferred protons (positive number 
for the transfer from 208 Pb to 40 Ca) is indicated as (a;p) 
(—6 < x < +5). The measured cross sections have been 
reported in Ref. [rJ 



cross sections towards the direction of the charge equilib- 
rium dominate. 

In (Op) and (-lp) panels of Fig. [HI the TDHF calcu- 
lation is seen to reproduce the measured cross sections up 
to six- neutron transfer. As the number of transferred pro- 
tons increases, (— 2p) to (— 6p), the cross sections in the 
TDHF calculation show a maximum at a neutron num- 
ber more than that of 40 Ca. Compared with measured 
cross sections, the TDHF results shift towards larger val- 
ues of neutron number. This behavior is similar to the 
case of 40 Ca+ 124 Sn reaction. Looking at the transfer 
cross sections for a fixed number of transferred neutrons 
in Fig. [T5j the TDHF calculations reproduce (+ln) and 
(+2n) panels rather well. 

As seen in Fig. [14l the TDHF calculations provide sub- 
stantial cross sections for proton pickup reactions, (+lp) 
to (+5p), which is the transfer towards the opposite di- 
rection of the charge equilibrium expected from the initial 
N/Z ratios. The cross sections show a peak around the 
neutron number 28. The TDHF calculations also provide 
substantial cross sections for many neutron pickup reac- 
tions (see bottom row of Fig. [T5|) . The cross sections show 
a peak around the proton number 20. These cross sec- 
tions come from an impact parameter region close to the 
fusion critical impact parameter. As seen in Fig. [121(a), a 
large average number of transferred neutrons up to 10 is 
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FIG. 15. (Color online) The same transfer cross sections for 
the 40 Ca+ 208 Pb reactions as those in Fig. 14. The number 
of transferred neutrons is indicated as (in) ( — 5 < x < +9). 



seen while the average number of transferred protons has 
small value. We note that the collision close to the fu- 
sion critical impact parameter accompanies large TKEL, 
and should suffer substantial evaporation effects which 
are not treated in the present analyses. 

The TDHF calculation systematically underestimates 
the cross section of neutron transfer processes from Ca 
to 208 Pb, (-In) to (-5n) (see top row of Fig. [T3J). Al- 
though these processes are against the charge equilib- 
rium, substantial cross sections are observed experimen- 
tally. In the TDHF calculation, cross sections of neutron 
transfer channels opposite to the charge equilibrium are 
several orders of magnitude smaller than the measure- 
ments. In Ref. jl9|, it has been argued that the neutron 
evaporation after collision is responsible for these chan- 
nels. 



C. 58 Ni+ 208 Pb reaction 

As a final case, we present results for the 58 Ni+ 208 Pb 
reaction at £ lab = 328.4 MeV (E c . m . ~ 256.8 MeV), for 
which measurements are reported in Ref. [llj ■ Since 
this system has ZpZy = 2296 exceeding the critical 
value 1600, we may expect an appearance of the quasi- 
fission process at a small impact parameter region. Us- 
ing the frozen-density approximation, the Coulomb bar- 
rier height is estimated to be Vb ~ 247.6 MeV, which is 
lower than the center-of-mass energy. We find the fusion 
critical impact parameter bf given by 1.38 fm for this re- 
action. To decide whether the nucleus once gets fused 
eventually decays into fragments or not, we continue to 
calculate the time evolution up to 4000 fm/c after two 



nuclei touches. If the fused system keeps compact form 
for this period, we regard the process as fusion. 



1. Overview of the reaction 

We first present an overview of the reaction dynam- 
ics. In Fig. [TH1 we show the deflection function in (a) 
and the TKEL in (b), as functions of impact parameter. 
In (a), we also show a deflection function for the pure 
Coulomb trajectory by a dotted line. The Coulomb rain- 
bow occurs at the impact parameter of 2.5 fm and the 
rainbow angle is 6* r — 121°. In Fig. [TTl we compare the 
Coulomb rainbow angle with measured differential cross 
sections which have been reported in Ref. [17[ . Red filled 
circles denote measured cross sections and blue solid ver- 
tical lines denote the Coulomb rainbow angle. From the 
figure, we find the measured differential cross sections 
show rather flat distributions compared with lighter sys- 
tems. This may be related to a rather small curvature 
of the deflection function around the Coulomb rainbow 
angle obtained from the TDHF trajectories as seen in 
Fig. HI (a). 

The TKEL shows a behavior different from lighter sys- 
tems. The maximum TKEL is about 50-60 MeV, similar 
to the value observed in 40 Ca+ 208 Pb reaction. However, 
there is a large impact parameter region, from 1.39 fm 
to 2.75 fm, in which the TKEL takes approximately the 
same value. 



58 Ni+ 208 pb (£ iab= 3 28 _4 MeV ) 
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FIG. 16. (Color online) Deflection function (a) and total ki- 
netic energy loss (b) as functions of impact parameter b for 
the reactions of 58 Ni+ 208 Pb at £i ab = 328.4 MeV. In (a), we 
show a deflection function for the pure Coulomb trajectory 
by a dotted line. 
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FIG. 17. (Color online) Differential cross sections of repre- 
sentative transfer channels as functions of scattering angle 
in the center-of-mass frame for the 58 Ni+ 208 Pb reaction at 
Slab = 328.4 MeV. The Coulomb rainbow angle obtained from 
the TDHF trajectories is denoted by blue solid vertical lines, 
and is compared with measured differential cross sections, red 
filled circles, which have been reported in Ref. [13]. 



In Fig. 1181 we show snapshots of density distribution 
for the 58 Ni+ 208 Pb reaction at the impact parameter of 
1.39 fm, just outside the fusion critical impact parame- 
ter. In the course of the collision, colliding nuclei form a 
rather thick neck and excurse for a long period connected 
by the neck. We find the two nuclei are connected for a 
period as long as 3600 fm/c. For collisions in the impact 
parameter region where the TKEL takes values around 
50-60 MeV, we find a formation of a similar thick neck 
which persists rather long period. These reactions are 
considered to correspond to the quasi-fission. 

Figure [19] shows the average number of transferred nu- 
cleons in (a), the N/Z ratios of the PLF and the TLF 
in (b), the average number of nucleons emitted to the 
continuum in (c) , and the fluctuation of the transferred 
nucleons in (d), as functions of impact parameter. 

In Fig. Q1J] (a), the average number of transferred neu- 
trons is shown by red filled circles connected with solid 
lines, while the average number of transferred protons is 
shown by green open circles connected with dotted lines. 
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FIG. 18. (Color online) Snapshots of density distribution 
of the 58 Ni+ 208 Pb reaction at E l!ih = 328.4 MeV and b = 
1.39 fm, just outside the fusion critical impact parameter. 



Positive numbers indicate the increase of the projectile 
nucleons (transfer from 208 Pb to 58 Ni) and negative num- 
bers indicate the decrease (transfer from 58 Ni to 208 Pb). 
From the figure, we find the average number of trans- 
ferred protons shows a minimum at b — 2.75 fm. We 
note that this value coincides with the impact parame- 
ter inside which the TKEL becomes almost constant in 
Fig. [T7J] (b). A similar minimum was also seen in the 
40 Ca+ 208 Pb case, as shown in Fig. GJ] (a). Outside this 
impact parameter, nucleons are transferred towards the 
direction of the charge equilibrium expected from the ini- 
tial N/Z ratios. In the impact parameter region, 1.55 fm 
< b < 2.75 fm, the average number of transferred nucle- 
ons increases as the impact parameter decreases for both 
neutrons and protons. A similar behavior was also seen 
in 40 Ca+ 208 Pb reaction as in Fig. [12] (a). 

At the impact parameter region b < 1.85 fm, the av- 
erage number of transferred protons becomes positive, 
opposite to the direction of the charge equilibrium of the 
initial system. However, the nucleon transfer still pro- 
ceeds towards the charge equilibrium of both the PLF 
and the TLF after the collision. This is clearly seen 
in Fig. [19] (b) which shows the N/Z ratios of the PLF 
and the TLF after collision. The N/Z ratio of the PLF 
(TLF) is denoted by red filled (green open) circles con- 
nected with solid (dotted) lines. As seen from the figure, 
the N/Z ratios of both the PLF and the TLF become 
closer to the N/Z ratio of the total system, 1.42, which 
is represented by a horizontal dashed line. 

As mentioned in the case of 40 Ca+ 208 Pb reactions, 
the change in the average number of transferred protons 
across the impact parameter b ~ 3 fm is related to the 
formation of the neck. Outside 6^3 fm, the neck is not 
formed and two nuclei are separated even at the closest 
approach. In such case, nucleons are transferred towards 
the direction of the charge equilibrium expected from the 
initial N/Z ratios. Inside 6^3 fm, the neck is formed 
between two nuclei. Then the transfer of nucleons pro- 
ceeds in two steps. Before the formation of the neck, the 
transfer of nucleons proceeds towards the charge equi- 
librium of the initial system in the same way as that in 
b > 3 fm. After the formation of the neck, an exchange 
of a large number of nucleons occurs at the time of the 
breaking of the neck. Depending on the position of the 
neck breaking, the transfer of nucleons is expected in ei- 
ther directions, from the target to the projectile or the 
reverse. Since the neck is formed with both protons and 
neutrons, the nucleon transfer in the neck breaking pro- 
cess accompanies both protons and neutrons in the same 
direction. 

Looking at Fig. Q1J] (a) , we find the increase of the av- 
erage numbers of transferred nucleons of both neutrons 
and protons as the impact parameter decreases below 
b = 2.75 fm. This indicates that the neck is broken 
at the position close to the target. Both protons and 
neutrons in the neck region are absorbed by the projec- 
tile. This mechanism explains the reason why the num- 
ber of transferred protons increases as the impact pa- 
rameter decreases in Fig. fTS] (a). This transfer process 
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FIG. 19. (Color online) The 58 Ni+ 208 Pb reaction at E lab = 
328.4 MeV. (a): Average number of transferred nucleons from 
the target to the projectile, (b): Neutron-to-proton ratios, 
N/Z, of the PLF and the TLF after collision, (c): Average 
number of nucleons emitted to the continuum, (d): Fluctu- 
ation of transferred nucleon number. The horizontal axis is 
the impact parameter b. In (b), the equilibrium N/Z value 
of the total system, 1.42, is indicated by a horizontal dashed 
line. 



associated with the neck breaking was also seen in the 
40 < 48 Ca+ 124 Sn and the 40 Ca+ 208 Pb reactions. 

At very small impact parameter region, 1.40 fm < b < 
1.50 fm, the average number of transferred neutrons 
shows a large fluctuation. The average number of trans- 
ferred protons also shows the fluctuation, correlated with 
that of neutrons. These fluctuations occur by changes of 
the breaking point of the neck. When the neck is bro- 



ken close to the target, a large number of nucleons are 
transferred from the target to the projectile, while the 
neck is broken at a midpoint between the projectile and 
the target, the number of transferred nucleons becomes 
small. 

Figure Q1J] (c) shows the average number of nucleons 
emitted to the continuum during the time evolution. The 
average number of neutrons (protons) emitted to the con- 
tinuum is denoted by red filled (green open) circles con- 
nected with solid (dotted) lines. As in other systems, we 
calculate it by subtracting the average number of nucle- 
ons inside a sphere of 14 fm for the TLF and that inside 
a sphere of 10 fm for the PLF from the total number of 
nucleons, 266. Again the number is rather small, about 
0.12 at the maximum. 

Figure [19] (d) shows the fluctuation of the transferred 
nucleon number. The fluctuation of transferred neutron 
(proton) number is denoted by red filled (green open) 
circles connected with solid (dotted) lines. They show 
a different behavior across the impact parameter around 
3 fm, indicating a qualitative change of the dynamics. 
Outside this impact parameter where protons and neu- 
trons are transferred in different directions, the fluctu- 
ation of transferred neutron number is larger than that 
of protons. Inside this impact parameter, although the 
average number of transferred neutrons is much larger 
than that of protons, the fluctuation is almost the same. 
This indicates that although the average number of trans- 
ferred protons is small, there is a strong mixture of single- 
particle orbitals of protons due to the formation and 
breaking of the neck. 

2. Transfer probabilities 

We next show transfer probabilities of the 58 Ni+ 208 Pb 
reaction as functions of impact parameter, which are 
shown in Fig. 1201 The small impact parameter region 
(b < 1.38 fm) corresponding to the fusion reaction are 
shaded. The positive (negative) number of the trans- 
ferred nucleons represents the number of nucleons added 
to (removed from) the projectile. The upper (lower) 
panels show neutron (proton) transfer probabilities for 
each transfer channel. In the left panels, we show trans- 
fer probabilities, (On) to (+6n) for neutrons and (Op) to 
(— 4p) for protons. They correspond to the transfer pro- 
cesses towards the charge equilibrium of the initial sys- 
tem. In the right panels, we show transfer probabilities, 
(+7n) to (+13n) for neutrons and (+lp) to (+5p) for 
protons, which dominate in the small impact parameter 
region, b < 2.75 fm. Probabilities of neutron transfer 
from 58 Ni to 208 Pb are very small and are not shown. 

In contrast to the previous cases of 40 ' 48 Ca+ 124 Sn and 
40 Ca+ 208 Pb, probabilities of transfer processes involv- 
ing more than 6 neutrons are seen in rather wide impact 
parameter region, 1.39 fm < b < 2.75 fm, where the for- 
mation of the thick neck is observed. In Fig. [TTJ] (b), a 
large value of TKEL was also seen in the impact param- 
eter region of b < 3 fm, indicating the significance of the 
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FIG. 20. (Color online) Neutron (upper panels) and pro- 
ton (lower panels) transfer probabilities as functions of im- 
pact parameter b for the reactions of 58 Ni+ 208 Pb at E\„\> = 
328.4 MeV. The positive (negative) number of transferred nu- 
cleons represents the number of nucleons added to (removed 
from) the projectile. Note that horizontal scales are differ- 
ent between the left and the right panels. A shaded region 
at small impact parameter (b < 1.38 fm) corresponds to the 
fusion reactions. 



evaporation effects. 

As in previous cases, probabilities of the processes ac- 
companying small number of exchanged nucleons show 
large spatial tail. The transfer probabilities for channels 
towards the charge equilibrium are large in most cases. 
The zero-proton transfer probability (Op, red solid line) 
in the left-bottom panel decreases as the impact param- 
eter decreases, show minimum at b ~ 3 fm, and again in- 
creases at smaller impact parameter region. This behav- 
ior is consistent with the behavior of the average number 
of transferred protons seen in Fig. [19] (a) . Although neu- 
tron transfer probabilities to the direction opposite to the 
charge equilibrium of the initial system are vanishingly 
small, we find appreciable probabilities of proton transfer 
opposite to the charge equilibrium of the initial system, 
as is seen from the right-bottom panel of Fig. [201 This 
feature is again consistent with the behavior of the aver- 
age number of transferred protons shown in Fig. [TO] (a) . 



3. Transfer cross sections 

We show transfer cross sections in Fig. |2"T1 and Fig. |2"21 
Each panel of Fig. EH shows cross sections classified ac- 
cording to the change of the proton number of the PLF 
from 58 Ni, as functions of neutron number of the PLF. 
Each panels of Fig. [22] shows cross sections classified ac- 
cording to the change of the neutron number of the PLF 
from 58 Ni, as functions of proton number of the PLF. Red 
filled circles denote measured cross sections and red solid 
lines denote results of the TDHF calculations. Again, re- 
action cross sections with relatively large values, such as 



1000 
100 

f ! 
is 0.1 
b 0.01 
0.001 
0.0001 

1000 
100 

f t 

0.1 



58 Ni+ 208 pb (£^-328.4 MeV) | TDHF 



(+5p)- 



(+4p)i 



(+3p) 



(+2p)i 



28 32 364044 28 32 36 40 44 28 32 36 4044 28 32 36 40 44 



I'M' 

(+lp) 



> 0.01 r [ V 

0.001 if \ 

o.oooi ' ■ 1 




28 32 3640 44 28 32 36 40 44 28 32 36 40 44 28 32 36 40 44 





(-5p)ir 



(-6p)i 



28 32 364044 28 32 36 4044 28 32 364044 28 32 36 4044 
NEUTRON NUMBER of PLF 

FIG. 21. (Color online) Transfer cross sections for the 
58 Ni+ 208 pb reaction at £ lab _ 328.4 MeV. Red filled circles 
denote measured cross sections and red solid lines denote re- 
sults of the TDHF calculations. The number of transferred 
protons (positive number for the transfer from 208 Pb to 58 Ni) 
is indicated as (arp) (—6 < x < +5). The measured cross 
sections have been reported in Ref. [l7| . 



(Op) and (-lp) panels of Fig. I2T1 and (+ln) and (+2n) 
panels of Fig. Ell are described reasonably well by the 
TDHF calculation. 

In Fig. 1211 as the transferred proton number increases, 
the calculation underestimates the measured cross sec- 
tion. The peak position of the cross section shifts to- 
wards larger neutron number compared with the mea- 
surements. A similar behavior was also seen in other sys- 
tems. This discrepancy is considered to be partly origi- 
nated from neutron evaporation processes which we have 
not yet taken into account. 

In (Op) and (-lp) panels of Fig.HH the TDHF calcula- 
tions overestimate the cross section for channels accom- 
panying large number of transferred neutrons (neutron 
number of PLF more than 34). We also find abundant 
cross sections for (+lp) to (+5p) processes, opposite to 
the charge equilibrium direction for the initial system and 
accompanying a large number of transferred neutrons. 
They come from reactions at small impact parameter re- 
gion, b < 3 fm, in which the transfer of nucleons associ- 
ated with the neck breaking is appreciable. 

In Fig. [22] we find an underestimation of cross sections 
for negative neutron transfer (— xn) (neutron transfer 
from 58 Ni to 208 Pb). On the other hand, almost constant 
cross sections are obtained for positive neutron trans- 
fer (+xn) (from 208 Pb to 58 Ni), up to the transfer of 10 
neutrons. The underestimation of the negative neutron 
transfer channels may be explained by the evaporation ef- 
fects as discussed in Ref. [17|. The cross sections for the 
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FIG. 22. (Color online) The same transfer cross sections for 
the 58 Ni+ 208 Pb reaction as those in Fig. 21. The number of 
transferred neutrons is indicated as (am) (—4 < x < +10). 



positive neutron transfer channels originate from reac- 
tions at small impact parameter which accompany large 
TKEL. Therefore, they may also suffer the evaporation 
effects. 



IV. COMPARISON WITH OTHER 
CALCULATIONS 

In this section, we compare our results of the TDHF 
calculations with those by other theories. Multinucleon 
transfer cross sections have been extensively and suc- 
cessfully analyzed by direct reaction theories such as 
GRAZING [35j and CWKB [36] . In both theories, rel- 
ative motion of colliding nuclei is treated in the semi- 
classical approximation. The probabilities of the mult- 
inucleon transfer processes are treated with a statisti- 
cal assumption using single-particle transfer probabilities 
evaluated with the time-dependent perturbation theory. 
We compare our results with those of the GRAZING 
for 40 ' 48 Ca+ 124 Sn reactions which have been reported 
in Refs. [HGJ]. 

In Figs. [23] and [24] we show transfer cross sections 
for the reactions of 40 Ca+ 124 Sn and 48 Ca+ 124 Sn, respec- 
tively. Each panel of these figures show cross sections for 
transfer channels classified according to the change of the 
proton number of the PLF as functions of neutron num- 
ber of the PLF. Red filled circles denote measured cross 
sections and red solid lines denote results of our TDHF 
calculations. Green crosses and blue open diamonds con- 
nected with dotted lines denote results of the GRAZING 
calculation. The latter symbols, blue open diamonds, 
show cross sections including the neutron evaporation ef- 
fect, while the former symbols, green crosses, without the 



evaporation effect. 

For the 40 Ca+ 124 Sn reaction shown in Fig. [23l one find 
that the cross sections by our calculation and those by 
the GRAZING are very close to each other for the pro- 
cesses shown in the panels of (Op), (— lp), (— 2p), (— 3p) 
and (— 4p). Cross sections accompanying many proton 
transfer, (— 5p) and (— 6p), are better described by the 
GRAZING compared with the TDHF. In both TDHF 
and GRAZING calculations, the peak positions of the 
cross sections are shifted towards large number of neu- 
trons in the (— 5p) and (— 6p) panels. The discrepancy is 
slightly remedied by including the neutron evaporation 
effect in the GRAZING calculation. 

For the 48 Ca+ 124 Sn reaction shown in Fig. [24] we again 
find a good coincidence between the TDHF results and 
those of the GRAZING. In both TDHF and GRAZING 
calculations, the cross sections shift towards the direction 
of small neutron number for (+2p) panel compared with 
measurements, while towards the direction of large neu- 
tron number for (— 2p) panel. In Ref. [ijj, the effect of 
the neutron evaporation has been evaluated to be small 
for this system. 

We notice that there are similar failures in the TDHF 
and the GRAZING calculations for the cross sections of 
channels accompanying transfer of large number of pro- 
tons. It seems that they are caused by a common prob- 
lem, although two theories are relied upon very different 
basis. One possible origin of the failure is an insufficient 
inclusion of the correlation effects beyond the mean-field 
theory. In the TDHF calculation, the many-body wave 
function is always assumed to be a single Slater deter- 
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FIG. 23. (Color online) Transfer cross sections for the 
40 Ca+ 124 Sn reaction at E^ h =170 MeV. Red filled circles 
denote measured cross sections, red solid lines denote results 
of the TDHF calculations, and green crosses (blue open dia- 
monds) connected with dotted lines denote calculated results 
using the GRAZING code without (with) the neutron evapo- 
ration effect. The number of transferred protons is indicated 
as {xp) (—6 < x < +1). The measured cross sections and the 
GRAZING results have been reported in Ref. [IoT |. 
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FIG. 24. (Color online) Transfer cross sections for the 
48 Ca+ 124 Sn reaction at Ei &h =174 MeV. Red filled circles 
denote measured cross sections, red solid lines denote results 
of the TDHF calculations, and green crosses (blue open dia- 
monds) connected with dotted lines denote calculated results 
using the GRAZING code without (with) the neutron evapo- 
ration effect. The number of transferred protons is indicated 
as (xp) (—2 < x < +2). The measured cross sections and the 
GRAZING results have been reported in Ref. Q]J. 



minant and correlations beyond the mean-field is not 
included. In the GRAZING calculation, multinucleon 
transfer probabilities are evaluated from single-nucleon 
transfer probabilities with a statistical assumption, ig- 
noring correlation effects among nucleons. 

We next consider an approach based on Langevin-type 
equations of motion which has been originally developed 
for and applied to fission dynamics [H, [89[ and has been 
recently extended to apply to multinucleon transfer reac- 
tions [3(| Hl|- We consider the 58 Ni+ 208 Pb reaction for 
which an application of the Langevin approach has been 
reported in Ref. [27$ . 

In the Langevin approach, multinucleon transfer pro- 
cesses are treated as sequential processes of single- 
nucleon transfers. In the theory, an empirical parameter 
describing nucleon transfer rate is introduced. Figure [55] 
shows cross sections for transfer channels of pure proton 
stripping without neutron transfer (left) and pure neu- 
tron pickup without proton transfer (right). Red filled 
circles denote measured cross sections, red solid lines de- 
note results of our TDHF calculations, and blue dotted 
lines denote results of the Langevin approach reported in 
Ref. [13. 

In the case of pure neutron pickup channels, (Op), the 
TDHF calculation gives a better description than the 
Langevin theory for cross sections up to four-neutron 
transfer. The TDHF calculation overestimates the cross 
sections for more than three neutrons, due to the quasi- 
fission process at small impact parameter region as dis- 
cussed in Sec. IIII CI On the other hand, cross sections of 
pure proton stripping channels, (On), are much better de- 
scribed by the Langevin theory than the TDHF, except 
for one-proton transfer channel. 

Figure [531 shows transfer cross sections for several pro- 
ton stripping channels. Again, red filled circles denote 
measured cross sections, red solid lines denote results of 
our TDHF calculations, and blue dotted lines denote re- 
sults of the Langevin approach. For these channels, the 
Langevin calculation gives a much better description for 
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FIG. 25. (Color online) Cross sections for transfer channels 
of pure proton stripping without neutron transfer (left) and 
pure neutron pickup without proton transfer (right) for the 
58 Ni+ 2os pb reaction at £ lab _ 328.4 MeV. Red filled circles 
denote measured cross sections [l7j . red solid lines denote 
results of the TDHF calculations, blue dotted lines denote 
results of the Langevin calculation [27| , and green open trian- 
gles connected with dotted lines denote results of the CWKB 
calculation 11711 - 



the transfer cross sections than the TDHF calculation. 
We should, however, note that an adjustable parame- 
ter describing the nucleon transfer rate is introduced in 
the Langevin approach, while no empirical parameter is 
introduced in the TDHF calculation once the Skyrme in- 
teraction is specified. In the calculation of the Langevin 
theory, evaporation effects are already included in the 
calculation. 

For the 58 Ni+ 208 Pb reaction, an analysis using the 
CWKB theory has also been reported in Ref. 



17J. In the 



CWKB theory, the multinucleon transfer processes are 
treated in a similar way to the GRAZING theory, evalu- 
ating statistically using single-nucleon transfer probabil- 
ities which are calculated by the first-order perturbation 
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FIG. 26. (Color online) Cross sections for transfer channels 
of (asp) (-3 < x < -1) for the 58 Ni+ 208 Pb reaction at E^b = 
328.4 MeV. The horizontal axis is the number of neutrons in 
the PLF. Red filled circles denote measured cross sections [13], 
red solid lines denote results of the TDHF calculations, blue 
dotted lines denote results of the Langevin calculation [2TT ] . 
and green open triangles connected with dotted lines denote 
results of the CWKB calculation ITtTI . 



20 



theory. In Figs. [25] and [26] the CWKB cross sections 
are shown by green open triangles connected with dotted 
lines. In Ref. [13], three results of cross sections have been 
reported: in a simple CWKB theory, adding the proton 
pair transfer effect, and taking account of evaporation ef- 
fects in addition to the proton pair transfer effect. We 
show in these two figures the simplest version of the cal- 
culation without the proton pair transfer effect and the 
evaporation. 

As seen from Fig.[3Sl the CWKB cross sections are very 
close to those of the TDHF except for AN > 5 in the (Op) 
panel. In Fig. [26] the CWKB cross sections are seen to be 
too small as the number of transferred protons increases. 
In Ref. [13] , proton pair transfer processes are introduced 
and added to the simple CWKB cross sections to examine 
the effect as a possible origin of the discrepancy. 



V. SUMMARY AND CONCLUSIONS 

In the present paper, we have reported fully micro- 
scopic calculations for the multinucleon transfer reac- 
tions in the time-dependent Hartree-Fock (TDHF) the- 
ory. We performed calculations for the reactions of 
40 < 48 Ca+ 124 Sn at E lah = 170, 174 MeV, 40 Ca+ 208 Pb 
at E lah = 235, 249 MeV, and 58 Ni+ 208 Pb at E lah = 
328.4 MeV, for which multinucleon transfer cross sections 
have been measured experimentally fiol [l3l li~7l Ii~9l | . We 
use a projection operator method [67[ to calculate the 
transfer probabilities as functions of impact parameter 
from the wave function after collision. From the reac- 
tion probabilities, we evaluate cross sections for various 
transfer channels. 

The systems we have calculated, 40,48 Ca+ 124 Sn, 
40 Ca+ 208 Pb, and 58 Ni+ 208 Pb show different behaviors 
in the multinucleon transfer processes characterized by 
the N/Z ratios of the projectile and the target, and by 
the product of the charge numbers, Z-pZ^. 

In the collisions with different N/Z ratios between the 
projectile and the target ( 40 Ca+ 124 Sn, 40 Ca+ 208 Pb, and 
58 Ni+ 208 Pb), we find a fast transfer of a few nucleons 
when the impact parameter is sufficiently large. The 
nucleons are transferred towards the direction of charge 
equilibrium expected from the N/Z ratios of the projec- 
tile and the target. This means that protons and neu- 
trons are transferred in the opposite directions. When 
the N/Z ratios are almost equal between the projectile 
and the target ( 48 Ca+ 124 Sn), we find a few nucleons are 
exchanged symmetrically. 

As the impact parameter decreases, a neck is formed 
at the contact of two nuclei. Then the transfer process 
proceeds in two steps. At the beginning of the reaction 
before the formation of the neck, a few nucleons are trans- 
ferred in the same way as described above. After forming 
the neck, the transfer of a number of nucleons occurs as 
a result of the neck breaking when two nuclei dissoci- 
ate. Since both protons and neutrons in the neck are 
transferred simultaneously in this mechanism, protons 
and neutrons are transferred in the same direction. 



As the charge number product ZpZ^ increases, there 
appears an impact parameter region where a thick neck 
is formed. We consider reactions in this impact parame- 
ter region correspond to the quasi-fission. As mentioned 
above, both protons and neutrons are transferred in the 
same direction when the neck is broken at the time of 
dissociation. A large energy transfer is also accompanied 
from the nucleus-nucleus relative motion to the internal 
excitations. 

Comparisons with measured cross sections show that 
the TDHF calculations describe cross sections reasonably 
well for transfer processes of a few nucleons between the 
projectile and the target. As the number of exchanged 
nucleons increases, the agreement becomes less accurate. 
When more than a few protons are transferred, cross sec- 
tions as functions of the transferred neutron number show 
a peak at the neutron number more than that in the 
measurements. The magnitude of the calculated cross 
sections becomes too small compared with the measure- 
ments. This discrepancy is expected to be, to some ex- 
tent, resolved when we introduce nucleon evaporation ef- 
fects in our calculations. 

We have compared transfer cross sections of the TDHF 
calculations with those by other theories. We find that 
results of the TDHF calculations are rather close to those 
of direct reaction model calculations such as GRAZING 
and Complex WKB. We should note that the Skyrme 
Hamiltonian used in the TDHF calculation is entirely 
determined from the ground state calculations and there 
is no parameter introduced to describe nuclear dynam- 
ics. We thus conclude that the fully microscopic TDHF 
theory may describe the multinucleon transfer cross sec- 
tions in the quality comparable to existing direct reaction 
theories. 

To increase the reliability of our calculation, an inclu- 
sion of evaporation effects is apparently important. Since 
the evaporation processes take place in much longer time 
scale than the reaction mechanism shown here, it is not 
realistic to achieve the calculation within the TDHF the- 
ory. Instead, we are considering to estimate them using 
a statistical model, putting excitation energy of the final 
fragments calculated in the TDHF theory as inputs. 

An inclusion of correlation effects beyond the mean- 
field approximation is also of much significance. At- 
tempts beyond the mean-field have been seriously un- 
dertaken recently in realistic calculations. For example, 
in Ref. [80| . calculations based on the Barian-Veneroni 
prescription have been reported. In Refs. [84l. l90l - [92l |. a 
stochastic effects are included in the mean-field dynam- 
ics. 

One of the most important correlations in nuclear low- 
energy dynamics is the pairing correlation. A full so- 
lution of the time-dependent Hartree-Fock-Bogoliubov 
(TDHFB) theory [93l - l96j would provide a satisfactory de- 
scription for it. At present, however, there have not yet 
been reported realistic TDHFB calculations of heavy-ion 
collisions. A simplified method to include the pairing 
effe cts i n the TDHF dynamics has also been discussed 
[97l4l01 |. Extending the present analysis including the 



21 



pairing correlation will certainly bring more satisfactory 
description and provide reliable understanding for the 
multinucleon transfer reactions. 
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Appendix A: Average and fluctuation of nucleon 
number 



We summarize formula for the average and the fluctu- 
ation of the nucleon number belonging to the PLF or to 
the TLF. Using the probability distribution P n , defined 



by Eq. (|6]) , the average and the fluctuation of the nucleon 
number belonging to the PLF are given by 



N 



n=0 



N 



\ n=0 



(Al) 



(A2) 



The average and the fluctuation of the nucleon num- 
ber belonging to the PLF may also be calculated directly 
from the orbitals without projection procedure. The 
average number is calculated from the density, p{r) = 



Y^i=i So- ^t( r i cr )V'i( r i integrating over the spatial 
region Vp: 



<M/p) 



dr p(r). 



(A3) 



The fluctuation can also be directly calculated from the 
single-particle orbitals by 



o = J ((Nv r ?) - (Nv,f 



\ 



N 



J/Vp 



(A4) 



[1] K. Sapotta, R. Bass, V. Hartmann, H. Noll, R.E. Ren- 

fordt, and K. Stelzer, Phys. Rev. C 31, 1297 (1985). 
[2] K.E. Rehm, A.M. van den Berg, J.J. Kolata, D.G. Ko- 

var, W. Kutschera, G. Rosner, G.S.F. Stephans, and 

J.L. Yntema, Phys. Rev. C 37, 2629 (1988). 
[3] R. Kiinkel, W. von Oertzen, B. Gebauer, H.G Bohlen, 

H.A. Bosser, B. Kohlmeyer, F. Puhlhofer, and D. Schiill, 

Phys. Lett. B 208, 355 (1988). 
[4] U.W. Scherer, W. Briichle, M. Briigger, C. Frink, H. 

Gaggeler, G. Herrmann, J.V. Kratz, K.J. Moody, M. 

Schadel, K. Summerer, N. Trautmann, and G. Wirth, 

Z. Phys. A 335, 421 (1990). 
[5] J. Speer, W. von Oertzen, D. Schiill, M. Wilpert, H.G. 

Bohlen, B. Gebauer, B. Kohlmeyer, and F. Puhlhofer, 

Phys. Lett. B 259, 422 (1991). 
[6] M. Wilpert, B. Gebauer, W. von Oertzen, Th. Wilpert, 

E. Stiliaris, and H.G. Bohlen, Phys. Rev C 44, 1081 

(1991). 

[7] W. von Oertzen, Z. Phys. A 342, 177 (1992). 

[8] GL. Jiang, K.E. Rehm, J. Gehring, B. Glagola, W. 

Kutschera, M. Rhein, and A.H. Wuosmaa, Phys. Lett. 

B 337, 59 (1994). 
[9] L. Corradi, A.M. Stefanini, D. Ackermann, S. Beghini, 

G. Montagnoli, C. Petrache, F. Scarlassara, C.H. Dasso, 

G. Pollarolo, and A. Winther, Phys. Rev. C 49, R2875 

(1994). 

[10] L. Corradi, J.H. He, D. Ackermann, A.M. Stefanini, A. 

Pisent, S. Beghini, G. Montagnoli, F. Scarlassara, GF. 

Segato, G. Pollarolo, C.H. Dasso, and A. Winther, Phys. 

Rev. C 54, 201 (1996). 
[11] G. Montagnoli, S. Beghini, F. Scarlassara, G.F. Segato, 



L. Corradi, C.J. Lin, and A.M. Stefanini, J. Phys. G 23, 
1431 (1997). 

[12] L. Corradi, A.M. Stefanini, J.H. He, C. Lin, S. Begh- 
ini, G. Montagnoli, F. Scarlassara, G.F. Segato, G. Pol- 
larolo, C.H. Dasso, and A. Winther, J. Phys. G 23, 1485 
(1997). 

[13] L. Corradi, A.M. Stefanini, J.H. He, S. Beghini, G. Mon- 
tagnoli, F. Scarlassara, G.F. Segato, G. Pollarolo, and 
C.H. Dasso, Phys. Rev. C 56, 938 (1997). 

[14] C.L. Jiang, K.E. Rehm, H. Esbensen, D.J. Blumenthal, 
B. Crowell, J. Gehring, B. Glagola, J. P. Schiffer, and 
A.H. Wuosmaa, Phys. Rev. C 57, 2393 (1998). 

[15] L. Corradi, A.M. Stefanini, C.J. Lin, S. Beghini, 
G. Montagnoli, F. Scarlassara, G. Pollarolo, and A. 
Winther, Phys. Rev. C 59, 261 (1999). 

[16] L. Corradi, A.M. Vinodkumar, A.M. Stefanini, D. Ack- 
ermann, M. Trotta, S. Beghini, G. Montagnoli, F. Scar- 
lassara, G. Pollarolo, F. Cerutti, and A. Winther, Phys. 
Rev. C 63, 021601(R) (2001). 

[17] L. Corradi, A.M. Vinodkumar, A.M. Stefanini, E. 
Fioretto, G. Prete, S. Beghini, G. Montagnoli, F. Scar- 
lassara, G. Pollarolo, F. Cerutti, and A. Winther, Phys. 
Rev. C 66, 024606 (2002). 

[18] S. Szilner, L. Corradi, F. Haas, G. Pollarolo, S. Begh- 
ini, BR. Behera, E. Caurier, E. Fioretto, A. Gadea, 
A. Latina, G. Montagnoli, F. Nowacki, F. Scarlas- 
sara, A.M. Stefanini, M. Trotta, A.M. Vinodkumar, and 
Y.W. Wu, Eur. Phys. J. A 21, 87 (2004). 

[19] S. Szilner, L. Corradi, G. Pollarolo, S. Beghini, BR. 
Behera, E. Fioretto, A. Gadea, F. Haas, A. Latina, G. 
Montagnoli, F. Scarlassara, A.M. Stefanini, M. Trotta, 



22 



A.M. Vinodkumar, and Y. Wu, Phys. Rev. C 71, 044610 [44 
(2005). 

[20] S. Szilner, C.A. Ur, L. Corradi, N. Marginean, G. Pol- [45 
larolo, A.M. Stefanini, S. Beghini, B.R. Behera, E. 
Fioretto, A. Gadea, B. Guiot, A. Latina, P. Mason, G. [46 
Montagnoli, F. Scarlassara, M. Trotta, G. de Angelis, [47 
F. Delia Vedova, E. Farnea, F. Haas, S. Lenzi, S. Lu- 
nardi, R. Marginean, R. Menegazzo, D.R. Napoli, M. 
Nespolo, I.V. Pokrovsky, F. Recchia, M. Romoli, M.-D. 
Salsac, N. Soic, and J.J. Valiente-Dobon, Phys. Rev. C [49 
76, 024604 (2007). 

[21] Sunil Kalkal, S. Mandal, N. Madhavan, A. Jhingan, [50; 
E. Prasad, Rohit Sandal, S. Nath, J. Gehlot, Ritika 
Garg, Gayatri Mohanto, Mansi Saxena, Savi Goyal, S. [51 
Verma, B.R. Behera, Suresh Kumar, U.D. Pramanik, 
A.K. Sinha, and R. Singh, Phys. Rev. C 83, 054607 [52 
(20011). 

[22] L. Corradi, S. Szilner, G. Pollarolo, G. Colo, P. Mason, [53; 
E. Farnea, E. Fioretto, A. Gadea, F. Haas, D. Jelavic- 
Malenica, N. Marginean, C. Michelagnoli, G. Montag- 
noli, D. Montanari, F. Scarlassara, N. Soic, A.M. Ste- [54 
fanini, C.A. Ur, and J.J. Valiente-Dobon, Phys. Rev. C 
84, 034603 (2011). [55 

[23] W. von Oertzen and A. Vitturi, Rep. Prog. Phys. 64, 

1247 (2001). [56' 

[24] L. Corradi, G. Pollarolo, and S. Szilner, J. Phys. G 36, 

113101 (2009). [57 

[25] C.H. Dasso, G. Pollarolo, and A. Winther, Phys. Rev. 

Lett. 73, 1907 (1994). [58 

[26] C.H. Dasso, G. Pollarolo, and A. Winther, Phys. Rev. 

C 52, 2264 (1995). [59 

[27] V. Zagrebaev and W. Greiner, Phys. Rev. Lett. 101, 

122701 (2008). [60 

[28] V.I. Zagrebaev and W. Greiner, Phys. Rev. C 83, 

044618 (2011). [61 

[29] S.C. Jeong, N. Imai, H. Ishiyama, Y. Hirayama, H. Miy- 

atake, and Y.X. Watanabe, KEK Rep. 2010-2 (2010). [62 

[30] V. Zagrebaev and W. Greiner, J. Phys. G 31, 825 

(2005). [63; 

[31] V.I. Zagrebaev, Yu.Ts. Oganessian, M.G. Itkis, and W. 

Greiner, Phys. Rev. C 73, 031602(R) (2006). [64 

[32] V. Zagrebaev and W. Greiner, J. Phys. G 34, 1 (2007). 

[33] V. Zagrebaev and W. Greiner, J. Phys. G 34, 2265 [65 
(2007). 

[34] C. Simenel, Eur. Phys. J. A 48, 152 (2012). [66 
[35] A. Winther, Nucl. Phys. A572, 191 (1994); A594, 203 

(1995). [67 
[36] E. Vigezzi and A. Winther, Ann. Phys. (N.Y.) 192, 432 [68 

(1989). 

[37] P.A.M. Dirac, Proc. Cambridge Philos. Soc. 26, 376 [69 
(1930). 

[38] P. Bonche, S. Koonin, and J.W. Negele, Phys. Rev. C [70 

13, 1226 (1976). [71 
[39] RY. Cusson, R.K. Smith, and J.A. Maruhn, Phys. Rev. 

Lett. 36, 1166 (1976). 
[40] S.E. Koonin, K.T.R. Davies, V. Maruhn-Rezwani, H. [72 

Feldmeier, S.J. Krieger, and J.W. Negele, Phys. Rev. C [73 

15, 1359 (1977). 

[41] K.T.R. Davies, V. Maruhn-Rezwani, S.E. Koonin, and [74 

J.W. Negele, Phys. Rev. Lett. 41, 632 (1978). 
[42] H. Flocard, S.E. Koonin, and M.S. Weiss, Phys. Rev. C [75 

17, 1682 (1978). 

[43] S.J. Krieger and K.T.R. Davies, Phys. Rev. C 18, 2567 [76 
(1978). [77 



K.T.R. Davies, K.R. Sandhya Devi, and M.R. Strayer, 
Phys. Rev. C 20, 1372 (1979). 

C.H. Dasso, T. D0ssing, and H.C. Pauli, Z. Phys. A 
289, 395 (1979). 

J.W. Negele, Rev. Mod. Phys. 54, 913 (1982). 

A.S. Umar, M.R. Strayer, and P.-G. Reinhard, Phys. 

Rev. Lett. 56, 2793 (1986). 

K.-H. Kim, T. Otsuka, and M. Tohyama, Phys. Rev. C 
50, R566 (1994). 

K.-H. Kim, T. Otsuka, and P. Bonche, J. Phys. G 23, 
1267 (1997). 

C. Simenel and Ph. Chomaz, Phys. Rev. C 68, 024302 
(2003). 

T. Nakatsukasa and K. Yabana, Phys. Rev. C 71, 
024301 (2005). 

A.S. Umar and V.E. Oberacker, Phys. Rev. C 71, 
034314 (2005). 

J.A. Maruhn, P.G. Reinhard, P.D. Stevenson, J. 
Rikovska Stone, and M.R. Strayer, Phys. Rev. C 71, 
064328 (2005). 

A.S. Umar and V.E. Oberacker, Phys. Rev. C 73, 
054607 (2006). 

J.A. Maruhn, P.-G. Reinhard, P.D. Stevenson, and M.R. 
Strayer, Phys. Rev. C 74, 027601 (2006). 
A.S. Umar and V.E. Oberacker, Phys. Rev. C 74, 
021601(R) (2006). 

Lu Guo, J.A. Maruhn, and P.-G. Reinhard, Phys. Rev. 
C 76, 014601 (2007). 

A.S. Umar and V.E. Oberacker, Phys. Rev. C 76, 
014614 (2007). 

Lu Guo, J.A. Maruhn, P.-G. Reinhard, and Y. 
Hashimoto, Phys. Rev. C 77, 041301(R) (2008). 
A.S. Umar, V.E. Oberacker, and J.A. Maruhn, Eur. 
Phys. J. A 37, 245 (2008). 

A.S. Umar and V.E. Oberacker, Phys. Rev. C 77, 
064605 (2008). 

K. Washiyama and D. Lacroix, Phys. Rev. C 78, 024610 
(2008). 

K. Washiyama, D. Lacroix, and S. Ayik, Phys. Rev. C 
79, 024609 (2009). 

C. Golabek and C. Simenel, Phys. Rev. Lett. 103, 
042701 (2009). 

D. J. Kedziora and C. Simenel, Phys. Rev. C 81, 044613 

(2010) . 

Y. Iwata, T. Otsuka, J.A. Maruhn, and N. Itagaki, 

Phys. Rev. Lett. 104, 252501 (2010). 

C. Simenel, Phys. Rev. Lett. 105, 192701 (2010). 

Y. Iwata and J.A. Maruhn, Phys. Rev. C 84, 014616 

(2011) . 

R. Keser, A.S. Umar, and V.E. Oberacker, Phys, Rev. 

C 85, 044606 (2012). 

Y. Iwata, J. Mod. Phys. 3, 476 (2012). 

C. Simenel, A. Wakhle, B. Avez, D.J. Hinde, R. du Ri- 

etz, M. Dasgupta, M. Evers, C.J. Lin, and D.H. Luong, 

EPJ Web of Conferences 38, 09001 (2012). 

T.H.R. Skyrme, Philos. Mag. 1, 1043 (1956). 

A. Bonaccorso, D.M. Brink, and L.Lo Monaco, J. Phys. 

G 13, 1407 (1987). 

A. Bonaccorso and D.M. Brink, Phys. Rev. C 38, 1776 
(1988). 

A. Bonaccorso and D.M. Brink, Phys. Rev. C 43, 299 
(1991). 

H. Freiesleben and J.V. Kratz, Phys. Rep. 106, 1 (1984). 
H. G&ggeler, T. Sikkeland, G. Wirth, W. Briichle, W. 



2.3 



Bogl, G. Franz, G. Herrmann, J.V. Kratz, M. Schadel, 
K. Siimmerer, and W. Weber, Z. Phys. A 316, 291 
(1984). 

[78] W.J. Swiatecki, Nucl. Phys. A376, 275 (1982). 
[79] R. Balian and M. Veneroni, Phys. Rev. Lett. 47, 1353 
(1981). 

[80] C. Simenel, Phys. Rev. Lett. 106, 112502 (2011). 
[81] H.J. Liidde and R.M. Dreizler, J. Phys. B 16, 3973 
(1983). 

[82] R. Nagano, K. Yabana, T. Tazawa, and Y. Abe, Phys. 

Rev. A 62, 062721 (2000). 
[83] J.W. Eastwood and D.R.K. Brownrigg, J. Compt. Phys. 

32, 24 (1979). 

[84] K. Washiyama, S. Ayik, and D. Lacroix, Phys. Rev. C 

80, 031602(R) (2009). 
[85] E. Chabanat, P. Bonche, P. Haensel, J. Meyer and R. 

Schaeffer, Nucl. Phys. A635, 231 (1998); A643, 441 

(1998). 

[86] K.A. Brueckner, J.R. Buchler, and M.M. Kelly, Phys. 

Rev. 173, 944 (1968). 
[87] V.Yu. Denisov and W. Norenberg, Eur. Phys. J. A 15, 

375 (2002). 

[88] Y. Abe, S. Ayik, P.-G. Reinhard, and E. Suraud, Phys. 



[90" 
[91 

[92 

[93" 
[94 

[95 

[96" 
[97 



[99 

[100 
[101 



Rep. 275, 49 (1996). 

P. Frobrich and I.I. Gontchar, Phys. Rep. 292, 131 
(1998). 

S. Ayik, Phys. Lett. B 658, 174 (2008). 

S. Ayik, K. Washiyama, and D. Lacroix, Phys. Rev. C 

79, 054606 (2009). 

B. Yilmaz, S. Ayik, D. Lacroix, and K. Washiyama, 

Phys. Rev, C 83, 064615 (2011). 

Y. Hashimoto and K. Nodeki, larXiv:0707.3083l [nucl-th]. 
B. Avez, C. Simenel, and Ph. Chomaz, Phys. Rev. C 
78, 044318 (2008). 

I. Stetcu, A. Bulgac, P. Magierski, and K.J. Roche, 

Phys. Rev. C 84, 051309(R) (2011). 

Y. Hashimono, Eur. Phys. J. A 48, 55 (2012). 

J. Blocki and H. Flocard, Nucl. Phys. A273, 45 (1976). 

S. Ebata, T. Nakatsukasa, T. Inakura, K. Yoshida, Y. 

Hashimoto, and K. Yabana, Phys. Rev. C 82, 034306 

(2010). 

G. Scamps, D. Lacroix, G.F. Bertsch, and K. 

Washiya ma, Phys. Rev. C 85, 034328 (2012). 

S. Ebata, larXiv:1211.6812l [nucl-th] . 

G. Scamps and D. Lacroix, Phys. Rev. C 87, 014605 

(2013). 



